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Abstract. Recent results on high transverse momentum (pt) hadron production in 
O | P+P, d+Au and Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) are 

■ reviewed. Comparison of the nuclear modification factors, RdAuiPr) and Raa(jPt), 

demonstrates that the large suppression in central Au+Au collisions is due to strong 
s ! ' final-state effects. Theoretical models which incorporate jet quenching via gluon 

Bremsstrahlung in the dense partonic medium that is expected in central Au+Au 
, collisions at ultra-relativistic energies are shown to reproduce the shape and magnitude 

C*~) ' of the observed suppression over the range of collision energies so far studied at RHIC. 
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There are now four years of beam data on high px hadron spectra available from the 
relativistic heavy ion collider (RHIC) covering multiple configurations (p+p, d+Au, 
Au+Au) and beam energies (62.4, 130, 200 GeV). The four RHIC experiments, 
BRAHMS, PHENIX, PHOBOS and STAR have observed and reported on several 
striking features from this new era of collider data, where the high cross-section for 
^ ■ pQCD hard-scattering processes plays an important role in heavy ion collisions for the 
first time. In particular, jets of high transverse momentum hadrons which come from 
the fragmentation of hard-scattered partons can be experimentally observed in heavy 
ion collisions using two-particle azimuthal correlations |2J- 

The yield of single inclusive high px charged hadrons leading these jets is observed 
to be substantially suppressed (by a factor of 3-5 for p T >6 GeV/c) in the most central 
Au+Au collisions at v / snn= 130 and 200 GeV compared to the expectation from scaling 
the yield in p+p collisions at the same beam energy by the average number of binary 
nucleon-nucleon collisions in Au+Au EJ EJ EJ- The yields in d+Au collisions are 
slightly enhanced at intermediate p? and consistent with binary-collision scaling for 

p T >7 GeV/c|Z| E| Ell- 
Detailed studies of the two-particle azimuthal correlations as a function of centrality 
in Au+Au collisions at a /snn= 200 GeV reveal that the correlations at dip = tc expected 
from the balancing partner jet created in the partonic hard-scattering become more 
and more suppressed in central Au+Au collisions in stark contrast to p+p and d+Au 
collisions. The back-to-back partner of the di-jet effectively "disappears" into the bulk 
matter soft hadronic background generated in the collision |21 E|. 
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Figure 1. Comparisons of h and it spectra from p+p(p) collisions to leading order 
and next-to-leading order pQCD parton model calculations. The left panel is from [T3] 
and the right panel is from |14j . 



These features, along with the observed large azimuthal anisotropy of hadron 
distributions at large pt in non-central Au+Au collisions'^], provide strong evidence 
that a dense, strongly interacting medium is formed in Au+Au collisions at RHIC, 
substantially effecting the production and dynamics of jets which traverse it. In this 
article I review the available data on the transverse momentum spectra of leading 
hadrons in p+p, d+Au and Au+Au collisions and how the theoretical interpretations 
of these data contribute to our understanding of the dense nuclear medium created in 
central heavy ion collisions at RHIC. 



2. The baseline: high p T particle spectra in elementary hadron collisions 

Within the perturbative QCD (pQCD) parton model, the short-distance perturbatively 
calculable hard-scattering cross-section factorizes from the long-distance non- 
perturbative parton distribution functions (PDFs) and parton fragmentation functions, 
providing a way to describe high p? particle production in elementary collisions|12j. The 
parton distribution functions are determined from deep inelastic scattering data, while 
the fragmentation functions are obtained from experimental studies of e+e- collisions. 
Within the pQCD parton model, the expression for the production of high transverse 
momentum hadrons in pp(p) collisions can be written 

j io = K 1^ dx a dx b f a {x aj Q )f b {x b ,Q ) , (1) 

d y d P T i&J nZc dt 
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where f a (x a , Q 2 ) and fb(xb, Q 2 ) represent the distribution of partons carrying momentum 

D° (z c ,Q 2 ) 

fraction x within the incoming hadrons, h/c — is the fragmentation function to 
convert parton c into colorless hadrons, h and do " (a ^- cd> is the hard-scattering cross- 
section, which is calculable in a perturbative series expansion in powers of the strong 
coupling constant, a s . Since the complexity of the calculation grows nonlinearly with 
the power of a s , leading order calculations with a phenomenological K-factor to account 
for the next-to-leading order corrections are often used. For v /snn>50 GeV, these 
calculations agree impressively well with the experimental data. In the left panel of 
Figure [TJ a parton model calculation is compared to data from pp(p) collisions at three 
collision energies [TB*]. The right panel of Figure ^ shows a comparison of both LO and 
NLO calculations to tt° production in p+p collisions at RHIC[14J. The addition of 
terms to account for finite transverse momentum Kt of the partons within the nucleus 
improves the agreement with the experimental data in both cases. The success of this 
description in describing the hadron spectra in p+p(p) collisions suggests that a similar 
formalism may be a useful starting point for the more complex systems: p(d)+A and 
A+A. 

3. Leading particle spectra in heavy ion collisions 

In p(d)+A and A+A collisions, additional nuclear medium effects must be taken 
into account. The parton distribution functions in nuclei differ from those in free 
nucleons[15j. The distributions of partons with x < 10 _1 are suppressed relative to 
the same distributions in nucleons, a phenomenon known as "shadowing" . For partons 
with x ~ 1 the distributions are enhanced or "anti-shadowed" . In addition, multiple 
soft scattering of the projectile partons as they traverse a target nucleus may boost 
their transverse momentum slightly even before they undergo the hard-scattering or 
subsequent to it, leading to an enhancement of the yield at high p T compared to p+p 
collisions [12]. These effects can be added to the factorized pQCD model expressed in 
Equation Q as modifications to the PDF terms. 

Finally, the nuclear medium produced in a collision might influence the production 
of high px hadrons, either through hard partonic or hadronic re-interactions in the final- 
state after the hard-scattering. Partonic interactions with a colored medium induce 
gluon Bremsstrahlung radiation, reducing the parton's energy before it fragments into 
hadrons [T7j. This effect is often introduced in pQCD parton model calculations through 
a modified fragmentation functional EJ] or in the form of quenching weights PU], for 
example. 

By comparing the yields in p(d)+A and A+A collisions to that in p+p collisions, 
with a scaling factor to take into account the nuclear geometry, one can test the 
assumption that a nucleus-nucleus collision is a simple superposition of incoherent 
nucleon-nucleon scatterings and explore these nuclear effects. The magnitude and shape 
of deviations from this scaling in Au+Au collisions may be used to infer information 
about the properties of the medium [TTH I21j. This is the fundamental goal of studying 
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Figure 2. Transverse momentum spectra and nuclear modification factor Raa of 
and 7r° measured by the PHENIX experiment for Au+Au collisions at v / snn= 130 
GeV0. 



high pt particle production and dynamics in heavy ion collisions: What can be learned 
about the medium from these "hard probes" ? 

The nuclear modification factor, Rab, is defined as 

R (rt h\ — d °AB/dyd 2 p T {b) ( . 

(binary (o) ) da£ p /dyd 2 p T 

where da h AB j 'dyd 2 pT{b) is the differential cross-section for hadron production in a nuclear 
collision of impact parameter, b, do-p p /dyd 2 pT is the corresponding quantity in p+p 
collisions and (iV binary (&)) is the average number of binary nucleon-nucleon collisions, 
often calculated using the Glauber model for the nuclear overlap with given impact 
parameter [22J. 

In the first collision data from RHIC at v / snn=130 GeV, the PHENIX experiment 
showed that the production of high pt hj^ and 7r°'s was suppressed by a factor of 
~3 in central Au+Au collisions compared to the expectation from binary-scaled p+p 
collisions, as seen in Figure El 0- This stunning result, also seen in spectra by the 
STAR experiment^), confirmed that a qualitatively new kind of system with intriguing 
properties is generated in these collisions. However these first data are quantitatively 
limited by the parameterized p+p reference spectrum used for the comparison. The 
reference had to be derived from a compilation of the world's data on hadron production 
in p+p(p) collisions over a broad range of collision energies because no data exist at 
./s^= 130 GeV. 

In RHIC Run II, this limitation was overcome when data from p+p and Au+Au 
collisions at v /snn= 200 GeV were collected |3 12H EDI- Figure 01 shows results for 
h Raa on the left and a similar measure, Rcp, the binary collision-scaled yield in 
central Au+Au compared to the scaled yield in peripheral Au+Au collisions [5] on the 
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Figure 3. Nuclear modification factors, Raa (left) and Rcp (right) for /i ± measured 
by the STAR experiment in Au+Au collisions at ^/snn= 200 GeVjS]. 



right. This comparison, which is motivated by the observed binary collision-scaling 
in peripheral Au+Au collisions for pt>2 GeV/c seen in the lower panels of the left 
side of Figure El extends the pt reach of the measurement with smaller systematic 
uncertainties than Raa- The high statistics data sample at v /snn= 200 GeV, with 
smaller uncertainties and higher pt reach than v /snn= 130 GeV, reveals the feature that 
for all hadrons with 6<pr<l2 GeV/c, the suppression is approximately independent of 
Pt- This pr-independence was largely unexpected in model predictions prior to the 
availability of the data [18]. 

More recent model calculations are compared to the data in Figure El Two 
formulations of a pQCD parton model incorporating partonic energy loss, labeled pQCD- 
I[25j and pQCD-II|2T] are shown. The influence of nuclear gluon shadowing is shown 
within pQCD-I as the dashed lines of the left panel of Figure El It is the subtle interplay 
of the shadowing, initial state multiple scattering and jet quenching which results in the 
observed p^-independence in these models. In addition, the right panel of Figure El 
also shows a calculation from the parton saturation model which extends initial-state 
gluon saturation beyond the saturation scale, Q s to high pr at midrapidityj2E]- This 
model is also seen to reproduce the suppression pattern for p T >5 GeV/c in the upper 
panel. The fact that two models based on entirely different physics, one incorporating 
strong initial-state effects and one with strong final-state effects, can both explain the 
suppression data was the primary reason that d+Au collisions were run at RHIC. Initial- 
state nuclear effects are present in both d+Au and Au+Au collisions, while final-state 
effects are expected only in Au+Au. 

Figure |H shows results on RdAu from each of the four RHIC experiments [Tj I%1 F)\ ITU] . 
The enhancement observed for 2<pt<Q GeV/c is similar to the enhancement seen 
in p+A collisions at lower energies f?Z], called the "Cronin effect". The standard 
explanation of this observed enhancement is the multiple soft interactions prior to and 
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Figure 4. Nuclear modification factors RdAu in d+Au collisions at • v /snn= 200 
GeV. Clockwise from the upper left: STAR h*®, PHENIX 7r°[B], PHOBOS ^[7] 
and BRAHMS ^(TU]. 
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Figure 5. A comparison of results from the Hadron-String-Dynamics modeljjSI] to 
RHIC data on Raa m Au+Au collisions at ^/snn= 200 GeV. 



possibly after the hard-scattering which were discussed in the first section [Tij]. 

In contrast, the gluon saturation model shown in the right panel of Figure El 



High pr hadron spectra at RHIC: an overview 



7 



predicted a suppression for central d+Au collisions of R^a — 0.7j2Hl- Clearly no such 
suppression is present in Figure EJ The factor 4-5 suppression in central Au+Au 
collisions must be attributed to final-state, not initial-state nuclear effects. 

Although the data in Figure H3 compare quite well with models incorporating 
final-state partonic energy loss, final state hadronic re-scattering may play some role 
in generating the observed suppression and many theoretical investigations of this 
possibility have been undertaken |281 I3T] . In Figured a comparison is made between 
RHIC data and a Hadron- String-Dynamics (HSD) model which incorporates quark, 
diquark, string and hadronic interactions with the medium 30J. The features of the 
data are indeed reproduced, but hadronic interactions are found to contribute very 
little to the suppression for p T >6 GeV/c, where inelastic, possibly colored partonic 
interactions of leading "pre-hadrons" with the dense medium are responsible for the 
observed suppression. 

4. Identified particle spectra and pseudorapidity dependence of Raa 

The particle composition at intermediate pt has been observed to contain an unusually 
bountiful baryon population in Au+Au collisions at RHICjSSESj- Figure El shows results 
from PHENIXlH] on the left and STAR[3lj on the right for Rcp of identified hadrons. 
The difference between meson and baryon Rcp is most pronounced for 2<p^<6 GeV/c. 
Several models relying on parton recombination and coalescence have been proposed to 
explain the "baryon anomaly" at intermediate From the right panel of Figure 

EJ it appears that the mesons and baryons converge for pt>Q GeV/c, where parton 
fragmentation dynamics dominate and no difference between mesons and baryons is 
expected. Interesting dynamics are clearly at play in the intermediate pt region, which 
should be investigated more thoroughly. What is also clear is the importance of carrying 
all leading hadron observables, such as spectra, di-hadron correlations and azimuthal 
anisotropies to the highest px available and at minimum for pt>Q GeV/c. This ensures 
that they can be used as reliable hard probes of the medium properties while minimizing 
complications from the dynamics of additional effects. 

In Figures [7| and |HJ data from the BRAHMS collaboration on Raa and Rcp 
at forward rapidities are compared to the results at mid-rapidity in Au+Au JU] and 
d+AujnS! collisions. From Figure [3 it appears that there is no strong evolution of the 
suppression pattern as a function of pseudorapidity in Au+Au, while in Figure |H1 there 
is a clear change from Cronin enhancement at midrapidity to a suppression at forward 
rapidities in the Au-going direction of d+Au collisions. Such a pattern is expected in 
saturation models |HZl IHH1 [HO] but may not be a unique signature of gluon saturation. 
In fact, energy-momentum conservation in the hadron fragmentation region alone should 
result in a reduction of RdAu of about 30% between 77 = and rj = 3.2 [3T]. In addition, 
conventional nuclear shadowing corrections are also expected to result in a roughly 
30% decrease of RdAu over the same range of pseudorapidity. The authors of 31] find 
that they can reproduce the data over all rapidities incorporating these effects into a 
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Figure 6. Nuclear modification factor Rcp for identified hadrons from PHENIXjH] 
on the left and STAR|33) on the right. 




Figure 7. Nuclear modification factors Raa and Rc p in Au+Au collisions at ^/snn= 
200 GeV for midrapidity (right) and forward rapidity (left) measured by the BRAHMS 
collaboration^!]] . 



hadronic absorption model without invoking substantial gluon saturation in the initial 
wave-function of the nucleus. More studies of the pseudorapidity dependence of Rab 
and other high pt observables such as di-hadron correlations are probably needed before 
a concrete understanding of the underlying dynamics emerges. 

5. Evolution of Raa with beam energy: Enhancement vs. Suppression? 

A Bjorken model estimate of the initial energy density achieved in central Au+Au 
collisions at RHIC yields values 30 times that in cold nuclear matter ^0] while at the 
SPS the corresponding value is approximately a factor of 20 [H]. These values are both 
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Figure 8. Nuclear modification factors RdAu (top) and i?J^p (bottom) for four 
rapidity intervals in d+Au collisions at a /snn= 200 GeV measured by the BRAHMS 
Collaborationist] . 
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Figure 9. Preliminary results from the STAR Experiment on Raa and Rcp in 
Au+Au collisions at y/s^m— 62.4 GeV|35]. 



well above the factor 3-5 expected for the onset of deconfinement and yet the yield of 
7r° measured at the SPS appears enhanced by a factor of 2-4 relative to binary collision 
scaling jl2j. Is there a transition from enhancement to suppression somewhere between 
17< v /s^<130 GeV? 

During RHIC Run IV, Au+Au collisions at v /snn= 62.4 GeV were run to search 
for such a transition. The choice of beam energy was motivated by the fact that it lies 
roughly halfway between the SPS and RHIC top energies on a logarithmic scale and 
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there exists a large body of reference data from p+p collisions at the ISR at this beam 
energy [4*3*]. As it turns out, two key findings underscoring the importance of a reliable 
reference dataset have emerged from this exploration. 

First, a recent re-analysis of the parameterization used for the SPS Raa reference 
reveals large variations in the extracted Raa depending on the choice of data and 
functional form used to fit them In fact, a parameterization using the best fit 
to all of the available data results in Raa for 7r° with 2<p^<4 GeV/c that is consistent 
(within large uncertainties) with binary collision scaling at v / snn=17 GeV. Rqp may 
even be slightly suppressed at the SPSjUJ- 

Second, it is challenging to find a reference spectrum within the ISR datasetspSJ 
that has systematic uncertainties at the same level of precision as the RHIC data. Even 
compiling a composite spectrum from all available publications results in minimum 
systematic uncertainties of 20-30%. Therefore, the significance of Raa measurements 
at this beam energy may remain small until a reference run at RHIC can be performed. 
Nevertheless, given these large uncertainties, preliminary results presented at the 
RHIC/AGS Users Meeting in May 2004 show that even at ^fs^= 62.4 GeV there 
is substantial suppression for pt>6 GeV/c[H3 EBl El EHj- Figure El shows preliminary 
results on Raa (left) and Rep (right) for h^ 1 measured by the STAR experiment at 
midrapidity (|?7| < 0.5) and slightly forward (0.9 < < 0.5) [45 . The intermediate pr 
region is less suppressed at v /snn=62.4 GeV than 200 GeV, but for pt>6 GeV/c, the 
suppression is comparable and may also be approximately p?- independent. 

This beam energy and pr-independence of Raa is intriguing. One possibility is that 
it may be due simply to the interplay of the evolving power-law shape of the underlying 
partonic spectrum and the medium-induced energy loss for media of varying size and 
density 49 . For central Au+Au at v /snn=200 GeV, calculations with a time-averaged 
transport coefficient of g=5-15 GeV 2 /fm for the medium and a pQCD parton spectrum 
from PYTHIA reproduce the observed Raa for pr>5 GeV/c and are pr-independent 
even to the highest accessible pt- This has intriguing consequences for Raa at the LHC: 
the strong suppression in the denser, more opaque system expected in Pb+Pb collisions 
at v /snn= 5.5 TeV is balanced by the shallower partonic spectrum, which results in 
more high p? hadrons from hard-scatterings near the surface. The overall effect is Raa 
with similar shape and magnitude at RHIC and LHC, in contrast to other models [2"T|. 
It is suggested that Raa may only be able to set a lower limit for the system density 
when q > 4 GeV 2 /fm|3S]. 

At v /snn= 62.4 GeV, the transport coefficient is probably reduced somewhat and 
may offer more sensitivity to the medium properties, but again it competes with the 
steeper partonic spectrum. Figure [TU] shows Raa expected for v / snn= 62.4 GeV Au+Au 
assuming q ~ 3-5 GeV 2 /fmjini- The suppression is substantial, pT-independent and 
reasonably consistent with the preliminary data for pt>6 GeV/c and with several other 
recently published predictions [5 0). Clearly, the body of data collected at RHIC over 
many systems and beam energies provides important feedback for testing and refining 
these models so that reliable estimates of the medium properties may be extracted from 




Figure 10. The px dependence of Raa at high pr compared to RHIC data at ^/snn= 
200 GeV (left) and predictions for Raa at ^/snn= 62.4 GeV within the same model 
from Ref.gHl (right). 



them. 



6. Conclusion 



Systematic high-precision studies of single hadron inclusive spectra measured in p+p, 
d+Au and Au+Au collisions at RHIC have yielded several exciting observations: the px- 
independent factor 4-5 suppression for px>6 GeV/c, the complex composition of hadron 
species at intermediate 2<p T <6 GeV/c and the possible beam energy independence of 
Raa(.Pt>6GgV / c) for 62.4< v /snn<200 GeV. These features have stimulated extensive 
theoretical feedback to explain their origins. The body of data supports the idea that 
strong final-state partonic interactions are at work in the dense medium created in these 
collisions. In order to take this qualitative statement closer to quantitative evaluation 
of the medium properties, a comprehensive, coherent explanation of all the high p? 
phenomena: single inclusive hadron spectra, di-hadron correlations/modifications and 
azimuthal anisotropic emission at high px must be developed. As the simplest and most 
investigated of the observables, the single hadron spectra provide a strong baseline 
for these more detailed and differential explorations of the medium. The quality 
and quantity of the spectra will grow and improve as the RHIC program proceeds, 
transforming them from the pool of discovery they provided in the first four years of 
RHIC running to a tool of precision for future high px studies in heavy ion collisions at 
RHIC and beyond. 
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